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Abstract 

We summarise what lattice simulations have to say about the physical properties of con- 
tinuum SU(N) gauge theories in 3+1 dimensions. The quantities covered are: the glueball 
mass spectrum, the confining string tension, the temperature at which the theory becomes 
deconfined, the topological susceptibility, the value of the scale that governs the rate at 
which the coupling runs and the r parameter that characterises the static quark potential at 
intermediate distances. 



1 Introduction 



In recent years there has been a resurgence of interest in developing analytic techniques for 
calculating the physics of non-Abelian gauge theories. Examples include approaches using 
light cone quantisation ^ and the AdS calculations based on the conjecture by Maldacena 
0. In practice all such calculations make approximations and so, in order to evaluate the 
accuracy of the approach being used, comparison with known results is very useful. In this 
context 'known results' means whatever has been learned from the computer simulation of 
these theories. The purpose of this review is to provide a readily accessible summary for such 
theorists of what has been learned from lattice calculations about the continuum properties 
of non-Abelian gauge theories in 3+1 dimensions. 

The quantities I shall focus upon are the glueball mass spectrum, the confining string ten- 
sion a, the deconfining temperature T c , the A A j 5 parameter that sets the scale in perturbative 
calculations, a parameter called r$ that is a characteristic scale of the static quark potential 
at intermediate distances, and the topological susceptibility Xt- 

There exists substantial information on these quantities for SU(2) and SU(3) gauge theo- 
ries. For many of these quantities lattice calculations exist for a large enough range of lattice 
spacings that a controlled extrapolation to the continuum limit is possible. This continuum 
physics is summarised in Section 2. For some other quantities accurate calculations do exist 
but only for one or two values of the lattice spacing, a, and so a reliable continuum extrap- 
olation is not possible. Where such calculations have been performed for values of a that 
are very small, the lattice corrections to mass ratios should also be small. This applies in 
particular to some of the heavier states of the mass spectrum, and we list these in Section |3|. 
The reader should be cautious in using the latter numbers: most of them 'should' be close to 
their corresponding continuum values but they might not be. 

While SU(2) and SU(3) gauge theories are of obvious relevance, many analytic approaches 
find it simplest to calculate the N — > oo limit of SU(iV) gauge theories; and this limit is, in 
any case, of particular theoretical interest. So in the last section I will mention some very 
preliminary SU(4) calculations and discuss what all this tells us about the large- N limit. 

Although this review restricts itself to 3+1 dimensions, we note that from the point of 
view of most of these analytic approaches, non-Abelian gauge theories in 2+1 dimensions are 
equally interesting and challenging. Here there is a recent lattice calculation of the glueball 
mass spectrum and string tension, in units of the scale provided by the coupling g 2 , which 
provides accurate extrapolations to SU(iV = oo) There also exist calculations of T c for 
,D = 2 + 1 SU(2), SU(3) and SU(4) gauge theories §. 

For the sake of clarity I have relegated to appropriate Appendices the details of what I 
have done, and what lattice calculations I have used, in order to obtain the final continuum 
numbers that I quote in Section 2. I have tried to perform the continuum extrapolations in a 
coherent way throughout and I apologise to those authors who may feel that I have not done 
their calculations full justice in the process. For example, calculations with improved actions 
on coarse lattices raise particular difficulties; and, whenever in doubt, I have chosen to err on 
the conservative side. 



1 



2 SU(2) and SU(3) in the continuum 



The spectrum of SU(iV) gauge theories is composed of states that are composed entirely of 
gluons and so it is usual to call the single particle states (including resonances) by the generic 
name 'glueballs'. Their masses we will denote by mg or by mjpc where J PC denotes the 
quantum numbers of the state. These theories appear to be confining at low temperature, with 
a linear potential at large distances between static sources in the fundamental representation. 
We label the corresponding string tension by a. While this characterises the large-r part of 
the potential, one can characterise its behaviour at intermediate distances by the distance r 
at which the force, F, has a particular value. It has become customary to use the particular 
definition r 2 F(ro) = 1.65 (which corresponds to a value that can be calculated with precision 
on the lattice and which can be estimated with some reliability from the observed spectrum of 
heavy quark systems). At sufficiently high temperatures the theory must become deconfined, 
and we label the temperature at which this happens by T c . At short distances the theory 
becomes asymptotically free and the scale controlling the running of the coupling can be 
chosen to be Ajjg. SU(iV) gauge fields in 3+1 dimensions possess a topological charge, Q, 
and a correponding topological susceptibility \t = {Q 2 )/V in a space-time volume V. There 
are other condensates, in particular the gluon condensate, but we choose to stop here. 

In a lattice calculation the overall scale is set by the lattice spacing a which is determined 
implicitly by the choice made for the inverse bare coupling (3 = 2N/g 2 (a). So, for example, 
a mass calculation will produce a number for the quantity omcfa), where the first a tells us 
that the mass is in lattice units, and the other a indicates that the mass will be afflicted by 
lattice discretisation errors. We can rid ourselves of the former by calculating the ratio of two 
such masses: amc{a) / am' G (a) = mG{a)/m' G (a). For small a one can expand the lattice action 
in powers of a and hence show that the leading corrections are ||: 



where any convenient mass scale /z, such as mg itself, may be used and where the coefficients 
Ci . . . are power series in the (bare) coupling. 

In the results I present in this paper I have usually extrapolated lattice masses to the 
continuum limit using the string tension 



since that quantity is known with some accuracy over a very wide range of lattice spacings. I 
could have used r equally well and the reader can easily translate from one continuum ratio 
to another using the values I tabulate. (r /a is more accurately calculable than a^Jo and so 
is likely to become the scale of choice in the future.) 

I have attempted to include all the useful 'modern' lattice calculations of these physical 
quantities and have extrapolated them to the continuum limit using either a simple 0(a 2 ) 
correction, as in eqn(|2|), or using both a 2 and a 4 corrections. Details are in the Appendices. 
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m G /Vo- 


state 


SU(2) 


SU(3) 






q «4 _i_ n nq 
o.ifi zn u.uy 


o++* 




6.58 ± 0.36 


2++ 


5.62±0.26 


5.17 ±0.17 


2 ++* 




6.73 ±0.73 


o-+ 


6.53±0.56 


4.97 ±0.57 


2-+ 


7.46±0.50 


7.05 ±0.65 


1++ 


10.2±0.5 




1+- 




6.32 ±0.31 


1-+ 


[10.4±0.7] 




3++ 


9.0±0.7 




3-+ 


[9.8±1.4] 





Table 1: Glueball masses, in units of the string tension, in the continuum limit. Values in 
brackets have been obtained by extrapolating from only two lattice values and so should be 
treated with caution. 



m G /m ++ 


state 


SU(2) 


SU(3) 


0++* 




1.78 ±0.12 


2++ 


1.46 ±0.09 


1.42 ±0.06 


2 ++* 




1.85 ±0.20 


o-+ 


1.78 ±0.24 


1.34 ±0.18 


2~+ 


2.03 ±0.20 


1.94 ±0.20 


1++ 


2.75 ±0.15 




1+- 




1.73 ±0.09 


1-+ 


[3.03 ±0.31] 




3++ 


2.46 ±0.23 




3-+ 


[2.91 ± 0.47] 





Table 2: Glueball masses, in units of the lightest scalar glueball mass, in the continuum limit. 
Values in brackets have been obtained by extrapolating from only two lattice values and so 
should be treated with caution. 



3 





SU(2) 


SU(3) 




1.201 ±0.055 


1.195 ±0.010 




0.694 ±0.018 


0.640 ±0.015 






0.531 ±0.045 


xtN° 


0.486 ±0.010 


0.454 ±0.012 



Table 3: The scale r , the deconfining temperature, T c , the topological susceptibility, xi , 
and the (zero-flavour) A parameter in the MS scheme. All in the continuum and in units of 
the string tension. 



In Table |] I list the continuum glueball masses (in units of the string tension) that I 
obtain, in this way, for both SU(2) and SU(3). Another way to present this 'data' is as a 
ratio of glueball masses. This may be convenient for people who don't calculate the string 
tension. I present the masses in this way in Table 0. (For reasons explained in the Appendix 
these masses are not precisely what one would obtain by dividing the appropriate entries in 
Table |.) 

So much for what is known about the mass spectrum. In Table |3| I present the continuum 
values of the other promised quantities. 

The error estimates in these Tables are supposed to be realistic. They include, albeit often 
indirectly, not only the statistical errors but also some of the important systematic errors. 
For example, in the mass calculations sufficiently large lattices are used for the finite-volume 
corrections to be very small. Using large lattices is slower and so leads to larger statistical 
errors than using small lattices. Similarly the process of extrapolating to the continuum limit 
with one or two correction terms magnifies the statistical errors on the component lattice 
calculations. In this sense it would be misleading to say that the errors in the Table are 'only 
statistical'. The precautions taken during the calculation to control systematic errors typically 
enlarge the statistical errors so that one should think of the former as being included, albeit 
imperfectly, within the quoted errors. Perhaps the main caveat here concerns the identification 
of the value of the spin J. A given irreducible representation of the lattice rotational group 
will contain states with different values of J in the continuum limit. One usually labels the 
lattice masses by the minimum such value of J. For calculations using operators smeared 
over physical length scales this can usually be justified Q for the lightest states with given 
quantum numbers (most clearly for P = +). Thus nearly all the spin assignments in the 
Tables should be correct. Nonetheless this still leaves open the possibility that, for example, 
what we have called the ++ * is in reality the 4 ++ . In principle it is straightforward to deal 



with this ambiguity [[16] but, for now, the ambiguity remains. Apart from this caveat, it 
should be the case that the above numbers will withstand the test of time; at least within, 
say, 2 of the quoted standard deviations. 
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m G /Vo- 


state 


SU(2) ; a ~ 0.03/m 


SU(3) ; a ~ 0.05/m 


o ++ 


o n i n on 

3.60 ± 0.20 


o r o i n 1 o 

3.52 ± 0.12 


2 ++ 


5.55 ± 0.27 


5.16 ± 0.21 


i ++ 


c\ on _i_ n no 


n n _i_ n T 

9.0 ± 0.7 


nJ L 

3 ++ 


9.75 ± 0.85 


8.9 ±1.1 


n — h 
U 


f .U4 ± U.o ^ 


O.O ± U.O 


2-+ 


7.70 ± 0.50 


6.9 ±0.4 


1-+ 


10.3 ±0.9 




3-+ 


9.5 ± 1.7 




2+- 




7.7 ±1.0 


1+- 




6.6 ±0.6 


2~ 




8.9 ±0.8 


1— 




9.9 ± 1.1 



Table 4: Glueball masses in units of the string tension for fixed but small values of the lattice 
spacing, as described in the text. No SU(3) estimates are given for the ±_ , 3 ±_ , 1~ + , 3~ + 
glueballs which appear to be too heavy to give a useful signal. 



3 SU(2) and SU(3) at finite lattice spacing 

The most comprehensive SU(3) lattice spectrum calculation is to be found in 0. (This title 
will be usurped by [|S| once the latter ceases to be preliminary.) The calculations were on a 
32 3 64 lattice at (3 — 6.4, which corresponds to a lattice spacing a ~ 0.1216/-^/cr ~ 0.05/m. 
Similarly, in SU(2) there is a calculation M on a 48 3 56 lattice at /3 — 2.85 which corresponds 
to a lattice spacing a ~ 0.0636/^/cr ~ 0.03/m. The values in physical units come from using 
y/a ~ 450MeV ||10|| , and are introduced solely to make the point that these are very small 



lattice spacings; they should not be taken more seriously than that. 

We present the mass ratios from these papers in Table [|. Since the lattice spacings are very 
small it is reasonable to expect that the 0(a 2 ) lattice spacing corrections to the continuum 
values will also be small. However the reader should be aware that this represents an extra 
unquantified systematic error. It is clear that in those cases where the states also appear in 
Table [I], the reader should use the latter values. 

The reason that our knowledge of the glueball mass spectrum is relatively incomplete is 
that since the late 1980's there have been almost no dedicated glueball mass calculations; the 
few that have appeared have usually come as a side-product of other calculations. Fortunately 
this is now changing and we can hope to have accurate values for all the above states (and 
more) quite soon. For a fore-taste see ||, although these are still 'preliminary' and so have 
not been included in this review. 



5 



SU(od) 




3.56 ±0.18 


m 2 ++/V°" 


4.81 ±0.35 




0.597 ±0.030 




0.428 ±0.022 


r y/(J 


1.190 ±0.043 



Table 5: Various continuum mass ratios after extrapolation to N = oo from N = 2 and 
N = 3; assuming the validity of eqn((^) all the way down to N — 2. Errors correspond to 
X 2 = l- 




4 SU(4) and large N 

We see from Tables HJ- f| that the physical properties of SU(2) and SU(3) gauge theories are 
very similar (apart from the absence of any C = — states in the former case). It is tempting 
to see this as evidence that we are already close to the N = oo limit even with N = 2, 3. Since 
the leading corrections are expected to be 0{1/N 2 ) this suggestion is not absurd; and indeed 
it is exactly what is found to occur for SU(iV) theories in 2+1 dimensions |4]]. If we assume 
that this is the case then we can extrapolate mass ratios to the N = oo limit using 

(3) 

and then we obtain the values displayed in Table ||. 

Of course, since the fits have no degrees of freedom, we have no information on the reli- 
ability of using eqn(^|) all the way down to SU(2), and the numbers in Table [5] must remain 
speculative. Clearly we need calculations for at least one further value of N; in that case we 
could test the goodness of the fits based on eqn(0). Unfortunately, while some exploratory 
SU(4) calculations do exist (JTI]], as summarised in Table [], they are not accurate enough to 
permit a useful continuum extrapolation. However all is not lost: the large- N arguments 
apply not only to the continuum limit but also to the theory at finite values of the lattice 
spacing. This enables us to squeeze some interesting information from Table |] as we shall 
now see. 

The question we want to ask is whether there is a smooth large- N limit and whether this 
is obtained by varying the coupling as g 2 oc 1/N. This is what one expects from the usual 
analysis |T2| of Feynman diagrams (and what one finds to occur in 2+1 dimensions [|J). 

Since the coupling runs the latter part of the question needs some elaboration. The usual 
/3-function will ensure that g 2 oc 1/N for large iV provided that Aj^ becomes independent 
of iV when expressed in units of a physical length scale in the SU(iV) theory, such as y/a. 
Equivalently suppose we define the running coupling, g%(l), on a length scale I = c\^ a where 
£ CT = \j\/o and q is small enough for the coupling to be small. Then what we are looking for 
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ST/(4) 




(3= 10.7 


P = 10.9 


/3 = 11.1 


a^Jo 
am ++ 
am ++* 
am 2 ++ 


0.296 ±0.015 
0.98 ±0.17 
1.54 ±0.44 
1.28 ±0.27 


0.228 ±0.007 
0.75 ±0.07 
1.39 ±0.16 
1.21 ±0.11 


0.197 ±0.008 
0.77 ±0.06 

1.03 ±0.14 

1.04 ±0.12 



Table 6: SU(4) masses and string tensions calculated in [11], on 10 4 , 12 4 and 16 4 lattices at 
(3 =10.7, 10.9 and 11.1 respectively. 



is 

^^) = ^„^) + 0(^) (4) 

for all N larger than some iVo- At the same time we want physical mass ratios, bar some 
obvious exceptions, to have finite non-zero limits as iV — > oo. 

We can test eqn(|) using the SU(4) string tensions in Table We note that the bare 
coupling can be used to define a running coupling on the scale of the cut-off a: g 2 (a) = 2N/ (3. 
However it is well-known that this is a 'bad' way to define a coupling and that one can do 



much better with a mean-field [13] or tadpole-improved WM coupling defined by 



9M) = ^i(H) (5) 

where (TrU p ) is the average plaquette. We can extract three SU(4) couplings in this way 
corresponding to the three values of (3 in Table || The scales on which these couplings 
are defined are given by the corresponding values of a\fa. We now go to the SU(2) and 
SU(3) lattice string tension values tabulated in the Appendices and find (by interpolation) 
the mean- field improved lattice couplings corresponding to these three values of a^/a. We wish 
to test the hypothesis that Ng 2 is becoming independent of N at large N, so we list these 
couplings, multiplied by N, in Table 0. We observe that these values are indeed consistent 
with being independent of N. This provides some evidence that eqn(^) is indeed valid all 
the way down to SU(2). Presumably the 0(1/N) corrections would become visible in more 
accurate calculations. 

Of course this analysis only makes sense if the theory has a smooth limit at large N. We 
can investigate this at each value of a separately. In Table ||] we plot the values of m ++/^/a 
for iV = 2,3,4 at each value of ay/a. In Table | we do the same for m 2 ++/v^"- We observe 
that these ratios are indeed consistent with having a smooth, confining large- N limit although 
the errors are too large for this to be a taxing test. 

We conclude that there is significant evidence that D=3±l SU(iV) gauge theories are 'close' 
to SU(oo) all the way down to SU(2). However it is also clear that very much better SU(4) 
calculations are needed (and indeed are in progress). 
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Ng'm 


q = a^/a 


SU(2) 


SU(3) 


SU(4) 


0.197 ±0.008 
0.228 ±0.007 
0.296 ±0.015 


4.966 ±0.050 
5.116 ±0.036 
5.434 ± 0.076 


4.938 ±0.045 
5.097 ±0.036 
5.397 ±0.063 


4.943 
5.149 
5.396 



Table 7: Bare couplings multiplied by N on three different length scales I = q£ ct , each scale 
being the same for the SU(2), SU(3) and SU(4) gauge theories when expressed in physical 
units (here the string tension). 



m ++/ v / o" 




SU(2) 


SU(3) 


SU(4) 


0.197 ±0.008 
0.228 ±0.007 
0.296 ±0.015 


3.67 ±0.10 
3.62 ±0.07 
3.49 ±0.03 


3.33 ±0.05 
3.22 ±0.04 
2.94 ±0.08 


3.91 ±0.33 
3.28 ±0.30 
3.31 ±0.57 



Table 8: The scalar glueball as a function of N, for lattice spacings that are independent of 
N, in units of the string tension. 



m 2 ++/y/cr 




SU(2) 


SU(3) 


SU(4) 


0.197 ±0.008 
0.228 ±0.007 
0.296 ±0.015 


5.59 ±0.09 
5.58 ±0.03 
5.54 ±0.14 


5.05 ±0.03 
5.02 ±0.04 
4.92 ±0.21 


5.23 ±0.61 
5.31 ±0.48 
4.30 ±0.90 



Table 9: The tensor glueball as a function of N, for lattice spacings that are independent of 
N, in units of the string tension. 
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5 Comments and Prospects 



Although it might seem, from Table p], that we already have a great deal of information about 
the mass spectra of SU(2) and SU(3) gauge theories, this appearance is deceptive. As soon 
as one tries to match these spectra with predictions of theoretical calculations (for example 
H|, p|) or models (for example one soon realises that it is important to have predictions 



for several excitations of the same quantum numbers as well as for spins which are different 
in the continuum but which appear in the same representation of the lattice symmetry group. 
For example, J = and J = 4. Thus one needs to develop techniques |L6] that will resolve 



such different continuum spins, and one needs to be able to calculate accurately heavier excited 
states. The latter will clearly entail developing techniques for dealing with resonances and 
continuum states. 

In order to take the continuum limit one needs a substantial range of lattice spacings, a, 
for which accurate calculations are available; extending from relatively coarse a to as fine a 
as one can go. At the coarser end, the masses of heavier excited states are so large in lattice 
units, arriQ, as to be virtually incalculable. Here the use of actions with different timelike and 
spacelike lattice spacings will be crucial: a s can be large while a t rriG remains calculably small 
|T7| , As we have seen in Appendix B, given the statistical accuracy of such calculations, it 
is important that the renormalisation of the asymmetry parameter be determined explicitly. 
Since part of the calculation must be performed for coarse a, there is an argument for using 
improved actions which have smaller lattice discretisation errors. If one does so then one 
might be tempted to perform the whole calculation with coarse lattice spacings. However we 
have seen that even if such actions reduce the leading 0(a 2 ) errors, it is dangerous to assume 
that these are so small as to be negligible compared to the sub-leading 0(a 4 ) errors. That 
is to say, both terms need to be included in the continuum extrapolation. In that case it is 
useful to have calculations at smaller values of a to accurately fix the coefficients of such fits. 

In addition to mass spectra, we have reviewed calculations of several other interesting 
physical quantities. Here too there are questions that need answering. Most pressing is 
a convincing demonstration that the effective string theory that governs the long-distance 
physics of the confining 'flux tube' is indeed in the universality class of the simple bosonic string 
(leading to the usual Luscher correction in eqn@ and the corresponding correction for the 
static potential). Whether this is in fact so depends on the effective degrees of freedom (central 
charge) of the effective string theory. There is numerical evidence for the usual assumption but 
it is not yet accurate enough to be really convincing. Also needed is a resolution of the ~ 4% 
discrepancy between the two groups of calculations of the SU(3) deconfining temperature, as 
discussed in Appendix B. 

A major gap at present is the lack of an accurate SU(4) calculation. While we have 
managed to squeeze a surprising amount out of the crude lattice calculations fll| that are 



currently available, a proper quantitative discussion of the large-iV limit in 3+1 dimensions 
must await a much better set of calculations. 

Acknowledgements: I have benefitted from discussions with many of the authors I quote. 
I am grateful to Rainer Sommer for providing me with the unpublished SU(2) r values. 
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Appendix A : SU(2) 



In this Appendix I describe in some detail how the SU(2) continuum numbers in the main 
body of the text were arrived at. 

(a) the mass spectrum and string tension 

The most extensive calculations are for the string tension and the lightest scalar and tensor 



glueballs. The values that I have used for my continuum extrapolations are listed in Table [10 



There are also some calculations for other glueball quantum numbers and these are listed in 
Table [11]. All these calculations use the standard Wilson plaquette action [|18|]. The sources 



for these calculations are as follows: (3 = 2.0 and (3 = 2.1 are from |T9|; f3 = 2.2 to (3 = 2.6 



are from | 20fl ; except for the 32 4 lattice at /3 — 2.5 which comes from [^1] as does the (3 = 2.7 



calculation; (3 = 2.85 is from [^]. There are some other string tension calculations that I 
have not used. Notably @ where one finds ay/a = 0.1871(10), 0.1871(32) on 16 4 , 32 4 lattices 
at (3 = 2.5; ay/a = 0.1207(3) on a 48 3 64 lattice at = 2.635; and ay/a = 0.0911(4) on a 32 4 
lattice at /3 = 2.74. (This last lattice is probably too small, in physical units, to be entirely 
trustworthy.) These values are consistent with the ones we list. 

Some technical asides. The values listed are effective masses from t = a to t = 2a for 
(3 < 2.4 for the string tension, for (3 < 2.5 for the J = 0, 2 glueballs, and for all (3 for the 
J = 1, 3 glueballs. For other (higher) (3 I list effective masses from t = 2a to t = 3a. When 
I need a single value of the string tension at (3 — 2.4 or at (3 = 2.5, I use ay/a = 0.2673(15) 
and a^/a = 0.186(3) respectively. The string tensions are all obtained from the masses of 
torelons. These are loops of fundamental flux closed through the periodic spatial boundary 
and hence neither contractible nor, in a confining phase, breakable. For a loop of length L in 
lattice units, one obtains a from the mass mi of the loop using 

2 7T 

ami — a> . (6) 

3L 

The (leading) correction in eqn(^) is a universal string 'Casimir energy', and the coefficient 
used in eqn(|]) assumes that the effective confining string belongs to the simplest bosonic 
string universality class [2~3|. There is some numerical evidence for this apparently reasonable 



assumption but it is not so precise as to be convincing. This correction typically contributes 
~ 4% to the values in Table [l(J 

To obtain the continuum limits listed in Table [I] I have taken the lattice values in the above 
Tables and used eqn@. If the fit is bad then the lowest-/? values are dropped till a good fit 
is obtained. I use conservative 25% confidence level (CL) limits for the errors; less where the 
best fit has CL < 70%. In the case of the + and of the 2 + one can fit all the way down to 
quite coarse lattice spacings ((3 = 2.0 and (3 = 2.2 respectively) and so there I also fit with 
an extra 0(a 4 ) correction to check that the resulting continuum extrapolations are consistent 
with the values I quote. I follow a similar procedure to obtain the mass ratios in Table |2|. 

In addition to the above calculations, which use the Wilson plaquette action, there are a 
number of calculations using other actions. These provide a useful universality check and I list 



in Table 12 my extrapolations to the continuum limit of these lattice calculations. The first 
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SU(2) 


n 

P 


lattice 


a^Jo 


am + 


am 2 + 


am - 


am 2 - 


2.0 


6 d 12 


0.69(5) 


l. 76(16) 








2.1 


6 3 12 


0.608(16) 


1.74(13) 








2.2 


8 3 10 


0.467(10) 


1.44(8) 


2.76(45) 






2.3 


10 4 


0.3687(22) 


-i - a 1/1/ f \ \ 

1.216(23) 


1.94(8) 






2.4 


1 A 4 

14 


0.2686(21) 


0.959(20) 


"1 /j O / O O \ 

1.487(33) 






2.4 


16 4 


0.2660(21) 


0.936(23) 


1.520(30) 






2.5 


16 4 


0.1834(26) 


0.660(22) 


1.028(27) 






2.5 


20 4 


0.1881(28) 


0.723(23) 


1.051(26) 


1.241(63) 


1.34(13) 


2.5 


32 4 


0.187(12) 


0.70(4) 


0.99(4) 


1.37(14) 


1.42(11) 


2.6 


20 4 


0.1326(30) 


0.51(3) 


0.73(4) 






2.7 


32 4 


0.1005(20) 


0.39(3) 


0.55(3) 


0.63(6) 


0.72(4) 


2.85 


48 3 56 


0.0636(12) 


0.229(12) 


0.353(16) 


0.430(35) 


0.49(3) 



Table 10: Various SU(2) lattice calculations of the string tension and lightest J = and 
J = 2 masses. Errors in brackets. 



SU(2) 


P 


lattice 


am 3 + 


arris- 


ami+ 


ami- 


2.5 


32 4 


1.88(43) 




1.70(20) 




2.7 


32 4 


0.90(5) 


1.32(11) 


1.02(3) 


1.11(4) 


2.85 


48 3 56 


0.588(36) 


0.71(5) 


0.626(27) 


0.677(21) 



Table 11: Various SU(2) lattice calculations of the lightest J = 1 and J = 3 masses. Errors 
in brackets. 



three calculations are from [25|. All three use an anisotropic lattice action with the ratio of 



timelike and spacelike lattice spacings chosen to be small: £ = at/o-s 1- This is a trick [17] 
that enables one to work with coarse a s and yet have masses, a^mc, that are not so heavy as to 
be incalculable. The first two actions are variants of tadpole-improvement |L4| while the third 
is an (anisotropic) Wilson action. The calculations take advantage of the improvement by 
being performed on quite coarse lattices. For this reason we have extrapolated these results to 
the continuum with a more complicated C\a 2 a + C2{a 2 <j) 2 correction term. A potential problem 
with these anisotropic lattice calculations is that they have been performed with values of £ 
that vary strongly with (3. This makes it tricky to motivate a simple uniform continuum 
extrapolation. The Manton [^6j and Symanzik tree-level improved action || calculations are 
from f2~0fl . These calculations suffer from having been performed at only two values of (3, 
but these are large enough that a simple 0(a 2 ) correction should suffice. The errors on the 
extrapolations correspond to x 2 — 1 an d are not really reliable. Despite these caveats, it is 
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clear that we have significant confirmation of universality from all these results. 



SU(2) - improved actions 


action 




■m 2 +/\/cr 




m 2 -/y/a 


tadpolel(a s 7^ a t ) 
tadpole2(a s 7^ a t ) 

Wilson(a s 7^ a t ) 
Symanzik(a s = a t ) 

Manton(a s = a t ) 


3.97 ±0.31 
4.4 ±0.5 
3.76 ±0.22 
3.93 ±0.21 
4.01 ±0.28 


5.0 ±1.1 

5.91 ±0.44 
5.67 ±0.31 
5.72 ±0.25 


6.77 ±0.53 
6.71 ±0.56 


6.78 ±0.80 
7.49 ± 0.60 



Table 12: SU(2) mass spectra, in the continuum limit, from calculations using a variety of 
'improved' lattice actions (see text). 



(b) the r parameter 

The SU(2) r parameter has been calculated for a range of lattice spacings in [27j. In 



Table [T^ I present the raw numbers |28[] and the corresponding string tensions, which are 
taken from Table pi] either directly or by interpolation. One can extrapolate to the continuum 
limit using either an a 2 a or an (a/ro) 2 correction. These give essentially identical best fits: 
for example 

r 0v /a = 1.201(55) - 0.56(2.48)aV. (7) 

with 80% confidence level. 



SU(2) 


P 


r /a 


a^Ja 


2.50 


6.39 ±0.09 


0.186 ±0.003 


2.55 


7.55 ±0.11 


0.157 ±0.003 


2.60 


8.82 ±0.08 


0.1326 ±0.0030 


2.70 


12.05 ±0.22 


0.1005 ±0.0020 


2.85 


19.08 ±0.88 


0.0636 ±0.0012 



Table 13: Values of the r parameter and string tension. 



(c) the deconfining temperature T c 

A periodic Euclidean lattice of size N^N T with N s ^> N T may be used to calculate the 
static properties of the gauge theory at a temperature T = l/aN T . In Table O I list the 
critical coupling f3 c corresponding to the SU(2) deconfining temperature for various iV T (see 
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and references therein) together with my interpolated values of the string tension. (These 
values have been extrapolated to N s = oo using a finite-size scaling analysis.) Thus 

Tc = Auc^y • (8) 

and we can extrapolate to the continuum limit with 

T 

— - = 0.694(18) - 0.28(33) (aT c ) 2 (9) 
or using a a 2 a correction. The best fit has an excellent 90% confidence level. 



SU(2) 


N T 




a^J o 


2 


1.8800 ±0.0030 




3 


2.1768 ±0.0030 


0.4950 ±0.0110 


4 


2.2986 ±0.0006 


0.3704 ±0.0023 


5 


2.3726 ± 0.0045 


0.2920 ± 0.0045 


6 


2.4265 ±0.0030 


0.2435 ± 0.0032 


8 


2.5115 ±0.0040 


0.1787 ±0.0040 


16 


2.7395 ±0.0100 


0.0891 ±0.0034 



Table 14: Values of the deconfining critical coupling and string tension. 



(d) the topological susceptibility xt 

If Q is the topological charge of a gauge field then the susceptibility is defined as Xt — 
(Q 2 ) /V where V is the space-time volume. So on an L 4 lattice, a 4 Xt — (Q 2 )/^. How to 
calculate the topological charge Q of a lattice gauge field is a problem that is now well- 
understood ||10|| . Here we will summarise several calculations and use them to extract a 
continuum value for Xt /V a - Our description of the methods used is inadequately brief and 
the reader needs to consult the literature for answers to a variety of natural questions. 

The oldest reliable way to calculate the topological charge Q of a lattice gauge field is 
to locally smoothen the gauge field, for example by minimising its action density, and then 
to calculate a lattice topological charge density Ql(%) — a 4 Q(x) at every site, and finally to 
sum that up over all lattice sites. This approach goes by the name of 'cooling' [[30], [L(J. I 
will call the cooled total charge Q c . The fields will contain some charges that are narrow and 
these will be particularly influenced by the presence of the lattice cut-off. So it is useful to 
define a charge Q c < from which such narrow charges have been excluded. The corresponding 
two susceptibilites, Xc and Xc', will differ at finite-a but should converge in the continuum 
limit. (Note that, for the sake of clarity, I drop the subscript t on x-) We shall use a cut-off 
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on the instanton size that corresponds to all charges with a peak height QL{%peak) — 1/lQir 2 
being excluded. The calculations I shall use are from |^T[ |32], [lj|. In these calculations the 
number of cooling sweeps is typically 15 or 20. Since the topological charge varies negligibly 
with cooling in this range, we can safely average such slightly different calculations. A more 

~. An 



recent 'improved' cooling method, whose charges we label by Qi c , is to be found in | 
approach using a renormalisation group inspired form of cooling has been developed in [34 
and the corresponding charges we label Qrg- (For the latter we shall use the values after '0 
smears'.) There also exist methods that do not smoothen the gauge fields. In particular one 
can, instead, form smooth smeared operators and apply them directly to the original gauge 
fields, correcting explicitly for the additive and multiplicative renormalisations, as is done in 
35]. We shall refer to this susceptibility as Xsm- (We use the values after 2 'smears'.) Finally 
one can use a geometric form of the total topological charge [3~2"|. To avoid problems with 
dislocations this is calculated on blocked lattice fields and averaged over all 16 blockings. This 
charge we call Qg- (As one takes a — > the number of blockings will, in principle, need to be 
increased.) All these methods will give different results at a fixed value of a, because they deal 
differently with instantons whose sizes are near the cut-off. However they should converge to 
the same continuum value. 



In Table [15] I list the 6 different susceptibilities. It is apparent that the variation between 



the different calculations decreases with increasing (3. In Table [16] I give the continuum extrap- 
olation of Xt^/V* 7 f° r eac h of these six calculations. (I have used an a 2 a lattice correction.) 
They are all pretty much consistent with each other and the final value I infer, and quote in 
Table [3], is an average value x\^ I \/ a — 0.486(10). 



SU(2) 






1/4 

axj 


1/4 


1/4 

aXic 


1/4 

o-Xrg 


"Asm 


2.20 


0.1481(53) 












2.30 


0.1420(21) 


0.1024(25) 


0.1211(30) 








2.40 


0.1163(13) 


0.0937(27) 


0.1026(25) 


0.1216(18) 


0.1299(21) 




2.44 












0.1027(19) 


2.50 


0.0844(11) 


0.0760(12) 


0.0800(18) 


0.0913(30) 


0.0957(21) 




2.5115 












0.0820(11) 


2.57 












0.0681(12 


2.60 


0.0627(13) 


0.0585(14) 


0.0610(16) 


0.0605(46) 


0.0674(15) 





Table 15: The topological susceptibility, using various methods as defined in the text. 



Appendix B : SU(3) 



In this Appendix I will describe the SU(3) lattice calculations which I use to obtain the 
continuum numbers in the main body of this paper. 
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SU(2) : xTN° 


n = c 


n = d 


n = G 


n = Ic 


n = RG 


n = sm 


0.486(12) 


0.468(15) 


0.480(18) 


0.501(45) 


0.528(21) 


0.480(23) 



Table 16: The continuum topological susceptibility, from the values in Table [T5|, in units of 
the string tension. 



(a) the mass spectrum and the string tension 



In Table |17| I present the lattice calculations of the string tension that I use. Where there 
is more than one reference listed, I have taken an average of the different values (checking, of 
course, that they are consistent with each other). For later use, values of the string tension 
have also been estimated by interpolation for a variety of intermediate (3 values at which there 
are no direct calculations. 



P 



5.50 
5.54 
5.60 
5.70 
5.80 
5.83 
5.85 
5.90 
5.93 
5.95 
6.00 
6.07 
6.10 
6.17 
6.20 
6.30 
6.40 
6.50 
6.57 



SU(3) 



0.5830 ±0 
0.5727 ±0 
0.5064 ±0 
0.3879 ±0 
0.3176 ±0 
0.2991 ±0 
0.2874 ±0 
0.2613 ±0 
0.2485 ±0 
0.2391 ±0 
0.2188 ±0 
0.1964 ±0 
0.1876 ±0 
0.1684 ±0. 
0.1608 ±0 
0.1398 ±0 
0.1216 ±0 
0.1068 ±0 
0.0975 ±0 



0130 
0052 
0028 
0039 
0016 
0009 
0007 
0028 
0012 
0014 
0008 
0010 
0012 
0012 
0010 
0012 
0011 
0009 
0012 



sources 



int 
int 



n 



int 
int 



37j, H, gp} P, ffi 
int 
int 
int 



38], m m, n, n 

int 



41, 44 



m 



int 



Table 17: The string tension with sources ('int' means obtained by interpolation). 



The calculations of the SU(3) mass spectrum fall into two groups which are sufficiently 
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different that they need separate treatment. The first and older category, involves calculations 
that use the standard plaquette action and are performed down to small values of the lattice 
spacing a. These allow accurate continuum extrapolations for the lightest glueballs. For the 
heavier glueballs they are much poorer. This is partly because many of the earlier calculations 
did not attempt to calculate them and, as a result, there is often too short a lever arm for 
a controlled continuum extrapolation. In addition heavy glueballs have larger values of arriQ 
and so the correlation function will often drop into the statistical noise before an effective 
mass plateau is really clear - except at the very smallest values of a, where the lattices have 
to be very large and hence computationally expensive. 

The second category of calculations is recent and uses improved lattice actions so that 
useful calculations can be performed for coarse lattice spacings. To avoid large values of amc 
and the consequent danger of having no useful signal, they use an anisotropic action where 
the space and time lattice spacings differ and satisfy a t <C a s . In that case one can have a 
small and easily measured value of a t ma, even when the glueball is heavy, while at the same 
time a s is large so that the spatial lattice can be small. I will now discuss these two kinds of 
calculations with a s = a t and a s ^ a t , in turn. 

• a s = a t calculations 

The lattice mass values that I use are listed in Tables [18|, [H]. This is an (almost) complete 
collection of 'modern' mass calculations where systematic errors should be negligible. The 
calculations extend to small values of a, as we can see by substituting ■sja ~ 0.5 fm in Table |TT . 
Thus one expects that over at least part of the range of a in which calculations exist, the leading 
0(a 2 ) corrections should dominate and we can perform extrapolations using 

%(«) rac(O) , 2 , 1fVl 
. . . = . . . + ca a. (10) 

We find that this functional form fits the ++ lattice values for (5 > 5.7 and the 2 ++ lattice 
values for (3 > 5.9. For the h and 2 h we include all the lattice values in the fits, but for 
the l +_ we find that we have to drop the j3 = 5.9 value. For the ++ * we have no degrees 
of freedom so we have no feedback about the validity of the fit in that case. The resulting 
continuum values are shown, a little later, in Table |22|. 

• a s ^ a t calculations 

This type of calculation is very recent and so the important questions have not all been 
resolved in the way that they have for the older and now well-established a s = a t calculations. 
On the other hand these calculations certainly embody features that will be part of future 
mass calculations, so it is worth discussing them in some detail. 

At the moment there is only one complete published glueball spectrum calculation of this 
kind ^7|. To be precise, this paper contains two separate calculations, one with a tree- level 
anisotropy of £ = {a s /a t )tree — 3 and a second with £ = (a s /a t ) tree = 5. An improved action 
is used which should reduce the 0(a 2 ) errors down to 0(a s a 2 ). The 0(a 2 ) errors should be 
small because a t is (very) small, and so are ignored. 

The lattice masses we use here are precisely the ones that are highlighted in the tables of 
p7| and there is no point reproducing all those readily accessible numbers here. However, to 
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5.50 
5.70 
5.70 
5.83 
5.90 
5.93 
6.00 
6.17 
6.20 
6.40 
6.40 



SU(3) : am G 



0+^ 



1.14 ± 
0.974 ± 

0.90 ± 
0.858 ± 
0.809 ± 
0.779 ± 
0.710 ± 
0.582 ± 
0.560 ± 
0.415 ± 
0.433 ± 



0.03 

0.029 

0.04 

0.043 

0.024 

0.008 

0.022 

0.010 

0.020 

0.014 

0.011 



2+" 



1.342 ± 0.039 
1.224 ±0.020 

1.11 ±0.03 
0.828 ±0.026 
0.845 ± 0.020 
0.609 ±0.016 
0.636 ±0.023 



0++* 



1.28 ±0.05 
1.00 ±0.03 



sources 



m 



35 1 

Hi 
Hi 



Hi 

381 



m 



E 



Table 18: The lightest 0+" 
lattice actions with a s = a t . 



and 2 ++ glueball masses and the first scalar excitation, using 



SU(3) : am G 


P 


o-+ 


2"+ 


1+- 


sources 


5.90 


1.65 ±0.20 


1.60 ±0.18 


1.90 ±0.20 




38 




6.00 


1.25 ±0.14 


1.28 ±0.15 


1.33 ±0.14 




38 




6.20 


0.92 ±0.08 


1.20 ±0.09 


0.984 ±0.045 




38 




6.40 


0.63 ±0.04 


0.81 ±0.04 


0.780 ± 0.030 




n 





Table 19: The lightest + , 2 + and 1 + glueball masses, using lattice actions with a s = a t . 



give the reader a flavour of their results we provide in Table ^U|the lightest ++ , 2 ++ (=E ++ ) 
and l + ~ masses together with the scale r in lattice units, all for the £ = 3 calculation. Also 
shown is the quantity r\a which, being a dimensionless ratio of physical quantities, is in a 
different category from the other quantities listed. Table |21| lists the same quantities for the 
£ = 5 calculation. We note that by using a t <C a s it is possible to get accurate lattice mass 
estimates even for the heavier states, despite the fact that the spatial lattice spacing is large: 
a s ~ 0.2 — 0.4/m. If we use the listed values of a s /ro and r^a and compare to the values 
in Table [17] we see that the smallest a s in these calculations corresponds to /3 ~ 5.60 in the 
symmetric lattice calculations. These are indeed very coarse (spatial) lattices. 

There are at least two questions that need answering if we are to be able to extract 
continuum numbers from these calculations: 

1) what is £ r = a s jap. 

2) what are the lattice corrections to the continuum limit of mass ratios? 
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SU(3) : C = 3 




a t m ++ 


a t m 2 ++ 


a t mi+- 




r 2 Q a 


1.7 


1.05(2) 


1.65(5) 


2.09(9) 


0.861(2) 


1.58(1) 


1.9 


0.87(1) 


1.47(4) 


1.80(6) 


0.773(2) 


1.52(2) 


2.0 


0.797(9) 


1.40(2) 


1.68(3) 


0.7271(8) 


1.462(7) 


2.2 


0.649(8) 


1.19(2) 


1.48(3) 


0.6192(8) 


1.362(8) 


2.4 


0.548(6) 


0.995(9) 


1.24(1) 


0.505(1) 


1.329(6) 


2.6 


0.464(7) 


0.781(8) 


0.97(1) 


0.4021(9) 


1.340(2) 



Table 20: The lightest ++ , 2 ++ (= E ++ ) and 1 + glueball masses, the parameter ro and its 
relation to the string tension, all for the £ = (a t /a s ) tr ee = 3 calculation. 



SU(3) : £ = 5 



p 


a t m ++ 


a t m 2 ++ 


a t mi+- 


a s /r 


r> 


1.7 


0.578(5) 


1.103(8) 


1.19(1) 


0.8169(9) 


1.473(9) 


1.9 


0.475(4) 


0.918(7) 


1.053(8) 


0.727(1) 


1.45(1) 


2.2 


0.362(3) 


0.686(4) 


0.819(4) 


0.5680(5) 


1.356(4) 


2.4 


0.303(3) 


0.542(2) 


0.652(5) 


0.459(1) 


1.342(4) 



Table 21: The lightest ++ , 2 ++ (= T 2 ++ ) and 1 + glueball masses, the parameter r and its 
relation to the string tension, all for the £ = {a t /a s )tree = 5 calculation. 



Let us start with the first question. The tree level value of a s /at = £ will be renormalised, 

'a s \ /a, 



at J \a t / 1 



x(l + 0(g% (11) 



tree 



The corrections here are potentially large and are not suppressed by powers of a, so they 
cannot be absorbed into the lattice corrections used in the continuum extrapolation of mass 
ratios. So, wherever £ r = a s /a t is important it should be calculated explicitly. This can be 
done straightforwardly (see Appendix D of [[| for example) if necessary. How necessary is it 
to do so? It is known ||8|, ^5j that £ r — £ f° r the type of improved actions used here. So one 
will have an idea of the spatial volume that is sufficient for control of finite volume effects, 
just from the tree-level estimate. Moreover one does not need to know £ r for ratios of glueball 
masses, since the scale a t disappears in such ratios and the correction terms will be affected 
in a minor way. However there is at least one place where it is important and that is for 
calculations of how the potential (which is obtained in units of at) varies with distance (which 
is in units of a s ). So we need to know £ r if we wish to express the glueball masses in units of 
r or ^Ja . The quantity r^a, on the other hand, is a dimensionless ratio of quantities derived 
from the potential and so the precise value of £ r should be less important for it. Since £ r is 
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not explicitly calculated in |57| this means that we have to be cautious about quantities that 
combine glueball masses with r or a. We shall see below that such caution is indeed justified. 

Now to the second question. The tadpole- improvement is supposed to ensure that the 
leading correction is 0(a s a 2 ) rather than 0(a 2 ). For calculations where a s is large, as it is 
here, one might argue that the notionally sub- leading 0(a*) correction will be much more 
important than the notionally leading 0(a s a 2 ) correction, so that the latter can be neglected 
and the continuum extrapolation can be performed using just the former. This has an obvious 
danger however: it will produce very small statistical errors at the price of potentially large 
systematic errors. 

There is no a priori obvious answer to this second question so we turn to what the cal- 
culations themselves tell us. The most accurately calculated mass is the lightest ++ and so 
we can ask what fit it prefers. We consider the dimensionless quantity r m ++ and try to 
perform a continuum extrapolation with corrections that are powers of a s /ro- For the £ = 3 
calculations we find that we can get a very good fit for (3 > 2.0, 



a, \ 2 /a s x 1 



r m ++ = 4.10(20) - (5.0 ± 1.2) ( — + (6.5 ± 1.6) — ) , (12) 



and poorer but still acceptable fits even if we include all the values in Table 20. For £ = 5 we 



can fit the values in Table |2l] equally well: 

r m ++ = 3.83(15) - (3.45 ± 0.75) ( — V + (4.50 ± 0.80) ( — V (13) 

Vr / \r J 

In neither case are statistically acceptable fits with just a 0(a 2 ) or just a 0(af) correction 
possible, as should be apparent from eqns (|T2] , [T3|) . Note also that the coefficient of the 0(a 2 ) 
term is not much smaller than that of the 0(af) term despite the fact that it presumably 
comes with an overall factor of a s . This is presumably because, for these very coarse lattice 
spacings, a s is not small. 

Although this already indicates that any continuum fit needs to incorporate both 0(a 2 ) 
and 0(af) correction terms, there is the possibility that the lightest ++ glueball is special. 
For example, its mass vanishes at the nearby critical point in the fundamental-adjoint coupling 
plane. (Although it is not apparent why this should lead to a larg 0(a 2 ) correction of the 
kind that we find.) Let us then consider the continuum extrapolation of r^a . For £ = 3 we 
obtain a good fit, for (3 > 1.9, as follows: 

r Q Ja = 1.192(10) - (0.32 ± 0.08) f—V + (0.66 ± 0.11) f—V (14) 

\r J \r J 

As we shall see below, this continuum limit is consistent with what one obtains from the 
standard a s = a t calculations. Once again fits with just a 0(a 2 ) or a 0(af) correction term 
are not possible. 

We infer from the above analysis that continuum extrapolations from this region of lattice 
spacings should include both 0(a 2 ) and O(a^) corrections. We can see from eqns ([l2| - |l4|) 



that the coefficient of the a 2 term is somewhat smaller than that of the at term, but that 
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this difference is not enough to make the a A s correction the dominant one. Thus we shall 
extrapolate all mass ratios using two correction terms as in eqn(|14T). 

In the above discussion we ignored the fact that we do not actually know the true value of 
a s /at. For the toa/cx extrapolation that does not matter much, but for r^m it does. The fact 
that it does is apparent from the fits in eqns([l2],[l3]). If we take those values and translate to 



the string tension using eqn(|T4D we find 

m ++ _ f 3.44(17) i = 3 
~ 1 3.21(13) £ = 5 



(15) 



Comparing with the value of 3.64 ± 0.09 in Table [3], it is clear that there is a problem with 
the values in eqn(|i~5|). The discrepancies are only ~ 5 — 10% but this is large at the level of 
statistical accuracy being achieved in these calculations. 

In conclusion, to minimise the systematic errors associated with the continuum extrapola- 
tion, we shall use two correction terms, as in eqn(|I~4"D, in all our continuum extrapolations. To 
minimise the error associated with deviations of £ r = a s / a t from the tree-level value, £, we shall 
extrapolate to the continuum limit either ratios of glueball masses, a t m G / a t m ++ = m G /m ++, 
or dimensionless ratios of quantities from the potential, such as r ^Jcr. To obtain the contin- 
uum value of mc/y/cr we multiply the extracted value of rriG/m Q ++ by the value of m ++/ v /a 
as obtained in the lattice calculations with a s = a t . If we do so we obtain the masses listed in 
Table [22|. (The a s ^ a t values are obtained by averaging the £ = 3,5 results, and for the 2 ++ 
we use both lattice representations.) Averaging these masses leads to the final mass estimates 
in Table [l] and Table [|. 



SU(3) 




m G /m ++ 




G 


a s ^ a t 


a s = a t 


a s ^ a t 


a s = at 


0++ 








3.64(9) 


2++ 


1.43(8) 


1.40(9) 


5.21(32) 


5.16(20) 


o-+ 




1.34(18) 




4.97(57) 


2-+ 




1.94(20) 




7.05(65) 


1+- 


1.70(10) 


1.81(17) 


6.19(39) 


6.56(52) 


0++* 


1.78(14) 


[1.77(21)] 


6.48(54) 


[6.65(47)] 


2 ++* 


1.85(20) 




6.73(73) 





Table 22: Continuum mass ratios for the two different kinds of calculation. Brackets denote 
continuum extrapolations from only two values of (3. 



To gain the full benefit of this kind of calculation, one clearly needs to calculate £ r =. a s /a t 
at every value of f3. This will have its own statistical error and is likely to increase the 
final errors on the masses by a significant amount. In addition it would be useful to have 
calculations down to smaller values of the lattice spacing (perhaps down to the equivalent of 
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(3 = 5.9 or (3 = 6.0 in the isotropic action calculations) so as to achieve a much better control 
of the Cid 2 + c 2 a A extrapolation. A calculation possessing these features would produce a very- 
accurate mass spectrum. 

(b) the r parameter 



In addition to the calculations of r^/a from |47| which were listed in Table ^Ojand Table [21 

uiette action |37|, f|9j and als 
O, H4J, M, |1. These are listed in Table || 



there have been two recent calculations using the Wilson plaquette action |37|, £|9| and also a 
scattering of older calculations, including 



SU(3) ; r /a 


P 


ref 4S] 


refQ37| 


ref 50, 40, 41, 44, 51| 


5.54 




2.054(13) 




5.60 




2.344(8) 




5.70 


2.922(9) 


2.990(24) 




5.80 


3.673(5) 






5.85 




4.103(12) 




5.95 


4.898(12) 






6.00 




5.369(9) 


5.35(3) 


6.07 


6.033(17) 






6.20 


7.380(26) 




7.37(3) 


6.40 


9.74(5) 




9.89(16) 


6.40 






9.75(17) 


6.50 






11.23(21) 


6.57 


12.38(7) 







Table 23: Values of r /a from calculations with the Wilson plaquette action. 



We now separately extrapolate to the continuum limit the lattice values in each of the three 
columns of Table [23], using an a 2 jr\ correction. (One gets almost identical fits using an a 2 a 
correction.) The continuum limits are listed in Table |24|. We also show the continuum limits 
for the anistropic calculations listed in Tables |20|, |21| . In this case we use a ci(a/r ) 2 + C2(a/r ) 4 
correction since, as we remarked earlier, we cannot get statistically acceptable fits with just 
an a 2 or an a 4 term. The £ = 5 calculation involves a fit with 3 parameters to 3 points and so 
we cannot tell how good a fit it is. I have therefore placed it in brackets. In the last column 
I show the range of (3 fitted in each well as the confidence level of the best fit and, in 

brackets, the CL corresponding to the error. 

We note that all the fits are very good apart from that using the calculations of ||37|| . One 
might imagine that this is because this calculation goes down to lower values of (3 where the 
0(a 2 ) correction becomes insufficient. However if we just fit the (3 > 5.70 part of the data we 
get a value of 1.179(13) with an even worse confidence level, ~ 20%. The fact that the fit is 
worse closer to the continuum limit is worrying and we therefore do not give this fit as much 



21 



weight as the three other fits. From these it seems that a final value of 1.195(10) is a safe one 
and so this is what we have quoted in Table [| 



SU(3) 




lim a ^ r \/a 


CL% 


range 


m 

[37] 

40, P@, 0, 51] 
[47| e = 3 
|47|£ = 5 


1.197(11) 

1.175(7) 

1.204(26) 

1.192(10) 

[1.177(13)] 


70(20) 
40(10) 
95(25) 
75(25) 

-(x 2 = i) 


5.70 <(3< 6.57 
5.54 < < 6.00 
6.0 < (3 < 6.5 
1.9 < < 2.6 
1.9 < < 2.4 



Table 24: Continuum values of r^a from various lattice calculations. 



(c) the deconfining temperature T c 

As described in Appendix A, a periodic Euclidean lattice of size N^N T with N s ^> N T may 
be used to calculate the properties of the gauge theory at a temperature T = l/aN T . Thus the 
deconfining temperature T c may be obtained at fixed N T by varying (3 and finding the value 
Pc{N T ) at which the system undergoes the deconfining phase transition. The corresponding 
temperature will be a{[3 c )T c = 1/N T in lattice units. There now exist calculations of this 
quantity using a variety of actions. 

The most extensive calculations are for the standard plaquette action and these are listed 
in Table 25. As well as the critical couplings, j3 c , I list my interpolated values of the string 
tension. The N T = 4, 6 calculations are from |52[ while the N T = 8, 12 calculations are from 
53fl . The values I give are those obtained after an extrapolation to infinite volume using the 



finite-size scaling relation [^, [51] 



/TV \ 6 

P c (N T ,N a )=p c (N T ,oa)-h(j±) . (16) 
For iV> = 4, 6 this relation has been fitted to a range of spatial volumes and the coefficient h 



explicitly determined [[5^, [54]]. For iV T = 8, 12 there is no such direct calculation but a value 
h = 0.06 ± 0.04 seems safe and this is what I have used in going from the values calculated 
in 



53fl for N a = 32, to the values listed in Table |25|. (Aside: the values of the string tension 
used in [^3] were, in retrospect, incorrect. This led to the continuum fit quoted there having 
an unacceptably low confidence level, and so it should be disregarded. This has also been 
noted in and their alternative continuum extrapolation is consistent with ours.) I have 
extrapolated these lattice values to the continuum limit using a (aT c ) 2 correction. The result, 
with confidence level, is given in Table 

There are in addition four other recent calculations that use different improved actions 
|55| , |56| , |57l |58| . I have extrapolated these separately to the continuum limit and the results 
are presented in Table [26|. I have used an 0(a 2 ) correction in all cases. If I use instead an 
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SU(3) 


N T 


/3 C (N T ,V = oo) 




4 


5.6925(2) 


0.395(4) 


6 


5.8941(5) 


0.265(2) 


8 


6.0618(11) 


0.1975(18) 


12 


6.3364(25) 


0.1325(13) 



Table 25: Values of the deconfining critical coupling, j3 c , and the string tension. 



0(a 4 ) correction, then the extrapolations do not move much - certainly not outside the errors 
quoted in Table 26. The main change is that the errors themselves are typically halved. 

What we see from Table ^ is that there is one group of calculations that favours a value 
T c /\/<7 ~ 0.630 and another group that favours T c /^/cr ~ 0.655. The statistical errors are 
small enough for these two values to be incompatible. It is not possible to make a completely 
rational choice between these two sets of calculations. However, given the much larger range 
of N T for the calculations with the plaquette action, and given the accuracy with which the 
string tension is known in that case, it is hard not view the set to which it belongs as being 
the more likely to be right. Nonetheless we wish to play safe, so we choose for Table ^ the 
value Tc/^/cr = 0.640 ± 0.015 which encompasses the various possibilities. 



SU(3) 


refs 


lim a ^ TJ y/a 


CL% 


range 


m 

\ 
\ 
\ 

\ 


52, 5^ 
55 
56 1 
57| 
58 1 


i 0.629(6) 
0.634(8) 
0.627(12) 
0.659(8) 
0.651(12) 


75(25) 
55(15) 

<x 2 = i) 

85(25) 
45(15) 


N T = 4,6,8, 12 
N T = 3,4,5,6 
N T = 4, 6 
N T = 2,3,4 
N T = 3,4,6 



Table 26: Continuum value of the deconfining temperature in units of the string tension, 
from the calculations as shown. 



(d) the topological susceptibility xt 

As for SU(2) there are a number of quite different ways of calculating the topological 
charge, Q, and the resulting topological susceptibility, Xt = (Q 2 )/V- The older calculations 
|59|1 use cooling; after typically 20 cooling sweeps they calculate the lattice topological charge 



density, Ql(x), whose sum over x provides an estimate for Q. Once again we can either 
take the full susceptibility, Xc or we can try to remove the narrow instantons that might be 
lattice artifacts. As in the case of SU(2) we shall present results where all instantons with 
Qiixpeak) > 1/167T 2 have been removed; the corresponding susceptibility we call Xd- (As for 
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SU(2) we suppress the subscript t on \ for clarity.) There is also a recent calculation that 
uses cooling |60[. As for SU(2) there are calculations that work directly with the Monte Carlo 



generated gauge fields, but use smeared operators and correct explicitly for the additive and 
multiplicative renormalisations [)61|. We call the resulting susceptibility Xsm- There are also 
calculations |62[] that use a smoothing of the gauge fields that is motivated by renormalisation 
group considerations. This gives Xrg- At /3 = 6.0 two of the above calculations provide values 
on two lattice sizes. The values of a\t^ corresponding to all these calculations are listed in 
Table E7L 



SU(3) 




axl /4 H 


axT H 


axi /4 m 


axl£ HI 


a X R G l 32 l 


5.70 


0.1595(23) 


0.1244(33) 








5.80 


0.1347(25) 


0.1147(28) 








5.85 










0.1179(28) 


5.90 


0.1135(17) 


0.1059(16) 




0.1019(21) 




6.00 


0.0915(21) 


0.0880(22) 


0.0977(34) 


0.0901(19) 


0.0937(12) 


6.00 






0.0870(40) 




0.0961(33) 


6.10 








0.0799(17) 


0.0816(20) 


6.20 


0.0722(36) 


0.0713(35) 


0.0701(25) 






6.40 






0.0537(48) 







Table 27: The topological susceptibility, using various methods as defined in the text. 



Extrapolating to the continuum limit using an 0(a 2 a) correction I obtain the continuum 



values shown in Table |28|. These are clearly consistent with each other. In addition, apart 
from the first two columns which come from the same configurations, these calculations are 
statistically independent (apart from the string tension). This provides the final continuum 
estimate in Table [3[ 



SU(3) : lim a ^o Xn/V a 




n = d @ 


n= c 


n = sm [|61|| 


n = RG[^ 


0.439(16) 


0.459(17) 


0.448(50) 


0.464(23) 


0.456(23) 



Table 28: The continuum topological susceptibility, extrapolating the values in Table f?7], in 
units of the string tension. 



(e) the A]yfg scale 

To calculate Aj^g one needs an accurate calculation of the running coupling over a large 
range of length scales. Some of the first realistic attempts involved calculations of the static 
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potential, with the running coupling being obtained from the Coulomb interaction at short 
distances. (See for example |pcf| .) The most recent calculation uses a different technique that 
allows a fine control of the various systematic errors and this is the value we have quoted 
in Table [3[ 
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